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Through an exact analysis using quantum Langevin dynamics, we demonstrate the crossover 
from ballistic to diffusive thermal transport in a harmonic chain with each site connected to Ohmic 
heat reservoirs. The temperatures of the two heat baths at the boundaries are specified from 
before whereas the temperatures of the interior heat reservoirs are determined self-consistently by 
demanding that in the steady state, on average, there is no heat current between any such (self- 
consistent) reservoir and the harmonic chain. Essence of our study is that the effective mean free 
path separating the ballistic regime of transport from the diffusive one emerges naturally. 
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Fourier's law is an old empirical law stating connection 
between heat current density and spatially varying tem- 
perature field. But it is still not clear what are the neces- 
sary and sufficient conditions for the validity of Fourier's 
law of heat transport [1, 2, 3]. Heat conduction through a 
one-dimensional ordered harmonic chain connected with 
two heat reservoirs at different temperatures shows bal- 
listic nature. Also it is well established that heat trans- 
port in one-dimensional momentum conserving systems 
(absence of external potentials) does not follow Fourier's 
law 0, § • Heat conduction in a long harmonic chain con- 
nected to self-consistent reservoirs at every site shows dif- 
fusive behaviour qualifying system size independent ther- 
mal conductivity @, 0, Q • The transition from ballistic 
to diffusive dynamics in thermal and electrical transport 
has recently received a lot of attention, fn a recent letter, 
Wang has reported to have obtained quantum ther- 
mal transport from classical molecular dynamics using 
a generalised Langevin equation of motion. Based on a 
"quasiclassical approximation" , the author claims to rec- 
oncile the quantum ballistic nature of thermal transport 
with diffusive one in a one-dimensional quartic on-site 
potential model. In Ref. [10( the authors have stud- 
ied the transition from diffusive to ballistic dynamics for 
a class of finite quantum models by an application of 
the time-convolutionless projection operator technique. 
Here through an exact analysis using quantum Langevin 
dynamics, we demonstrate the crossover from ballistic 
to diffusive thermal transport in a harmonic chain con- 
nected to self-consistent reservoirs. 

Consider heat conduction through a one-dimensional 
ordered harmonic chain of particles I — 1, 2. ..TV with unit 
masses which are connected by harmonic springs of equal 
strengths. The Hamiltonian of the system is 
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where {xi} are Heisenberg operators, correspond to parti- 
cle displacements about some equilibrium configuration. 
We choose the boundary conditions xq = xn+i — 0. All 
the particles are connected to Ohmic heat reservoirs with 



coupling strength controlled by dissipation constant 7;. 
We set ji = 7 for / = 1, N and 7; = 7' for I = 2, 3.JV- 1. 
This allows us to tune the coupling (7') between self- 
consistent reservoirs and the chain sites without affecting 
the couplings at the end reservoirs. The temperatures of 
the first and last reservoirs are fixed as T\ = Tr, and 
TV = Tr respectively. The temperatures of the attached 
interior reservoirs {T{\ are determined self-consistently 
by the condition of net zero heat current from the side 
reservoir to the chain. Slightly different version of this 
model has been studied in [7], l8[ for infinitely long chain 
length. The time-evolution of the chain particles is gov- 
erned by a combination of Hamiltonian and stochastic 
dynamics, where the nearest neighbour harmonic poten- 
tials form the Hamiltonic part and the stochastic influ- 
ence comes from the connection of each chain site with 
Ohmic heat bath. As pointed out in 0, the influences 
of the interior reservoirs can be considered as the effect 
of degrees of freedom not present in the Hamiltonian. It 
is also shown in [7[ that in the limit N — > 00, the steady 
state is a local equilibrium state with a temperature pro- 
file satisfying Fourier's law with a temperature indepen- 
dent, finite thermal conductivity for the classical model. 
The quantum version of this model is studied in Q where 
a finite, temperature dependent thermal conductivity is 
found in the quantum regime. 

The quantum Langevin equations of the chain sites are, 

x\ = -2xi + + xi + i - 7;x; + rji for I = 1, 2, 3..JV, 

where r\i is the noise generating from the /th reservoir. 
The correlations of noises, are such that the distribu- 
tions of normal modes in isolated reservoirs follow Bose- 
Einstein statistics. Noise-noise correlation in the fre- 
quency domain is given by, 

^(mi^Vmiv') + m(^')Vm{uj)) = ^ coth(^)6(w + u')6i 

with h = Kb = 1. Now our first task is to determine 
the temperature profile of the interior reservoirs from the 
self-consistent condition. For this we write down the heat 
current from the reservoir to the chain and then expand 




FIG. 1: Plot of the temperature profile {T;} as a function of FIG. 2: Plot of the temperature profile {Ti} as a function of 

scaled length l/N for different N with 7 = 1.0 and 7' = 0.1. scaled length l/N for different 7' with 7 = 1.0 and JV = 256. 

The inset shows temperature dependence of scaled current for The inset shows temperature dependence of scaled current for 

different N with above values of 7,7'. Here mean free path different 7' with above values of j,N. Here mean free path 

t = 30. t = 3/V. 



it in the linear response regime; the current from the ith 
reservoir to the chain is given as Q 
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where G(u>) is inverse of a tridiagonal matrix with offdi- 
agonal elements equal to —1 and diagonal elements equal 
to 2 — uj 2 — i^'uj except the ends which are equal to 
2 - uj 2 - iju. Also T = (T L + T R )/2. We set J; = 
for / = 2. 3...N — 1 and solve N — 2 linear equations nu- 
merically to find {Ti} profile. In FigQ]we plot {T/} for 
different lengths of the chain for some fixed small value of 
7'. In the limit 7' << 1 we find the temperature profile 
is scaled as 



T, 
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{T;} profile is linear. For larger value of 7', the temper- 
ature profile becomes linear for smaller system sizes. 

To find the heat current through the wire from the 
left to the right heat baths we can use the heat current 
expression of Eq.@ with I = 1 or N. But we notice that 
Eq.([2]) then requires the temperatures at the boundaries 
(2) accurately as these terms contribute significantly. So if 
we want to use the scaling form of the temperature which 
is not so good at the boundaries for smaller wire sizes, 
it is better to evalute current in the middle bond of the 
wire. Current through (I, I + 1) spring of the chain is 
given by Q 



(3) 



where I = 2>ji and AT = T B , - T L . Here 6 is the jump 
in the temperature at the boundaries. The above scaling 
relation can be derived from a persistent random walk 
model of phonons in analogy with one for electrons [111 ] ; 
here £ is interpreted as the mean free path of the phonons. 
We also plot {7]} profile in FigHfor fixed length N = 256 
with different values of 7'. When 7' tends to zero (then 
I goes to infinity) , the heat transport in the chain is bal- 
listic (becomes completely ballistic at 7' = 0), and the 
temperature profile is flat as shown in Fig[2j With in- 
creasing 7' system transits through a mixed transport 
regime towards a diffusive one for sufficiently large value 
of 7' where I is much smaller than the system size; then 
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Now using the numerical solution for {T}} we first evalute 
the heat current with varying temperature for different iV 
and 7', and plot it in the inset of Figs[TJand[2]respectively. 
•h,i+\ (call J) is independent of I and we calculate it in 
the bulk for accuracy. Using the scaling form of {TJ} we 
find 



J 



k(T) AT 
(N + i) 



(5) 



k(T) is the temperature dependent thermal conductivity 
of the infinite chain. Above current expression is exact 
for larger size of the chain but for smaller size there will 
be correction from the boundaries. It clearly shows that 
for N >> I transport is diffusive satisfying Fourier's law 
and in the opposite limit the current is independent of 
(ballistic) . We clarify that the cross-over from ballistic to 
diffusive behaviour in transport depends on the effective 



3 



length scale of the problem and can be controlled here by 
tuning £, i.e., 7'. This model has similarity to the one- 
dimensional quartic onsite potential model [9[| if one iden- 
tifies 7' with the strength of the quartic onsite potential. 
But in the quartic onsite potential model the temperature 
and the strength of quartic potential are conjugate to 
each other, i.e., for a fixed strength of the quartic poten- 
tial, increasing the temperature one can cross-over from 
ballistic to diffusive regime of transport; similarly for a 
constant temperature, changing the strength of quartic 
potential one can tune from ballistic to diffusive trans- 
port. But in case of the self-consistent reservoir model 
the temperature and the strength of the coupling to the 
reservoirs (7') are independent parameters, not affecting 
each other. 

In conclusion we have demonstrated both ballistic and 
diffusive regime of thermal transport within a single anal- 
ysis of quantum Langevin dynamics. This is contrary to 
the remark made by author in Ref . [9( . As discussed nicely 
in [H, @] , it is a big challenge to derive Fourier's law for 
a system from microscopic Hamiltonian bulk dynamics. 
One can think of the problem as either (a) the system 
in microcanonical ensemble evolving towards equilibrium 
from an initial arbitrary distribution and study the relax- 
ation mechanism from different correlations like energy- 
energy; or (b) the system is kept in a non-equlibrium 
steady state by connecting it at the boundaries with 
stochastic or mechanical reservoirs and then determine 
the size dependence of the steady state current . From 
a large number of numerical and a few analytical studies 
[l], it is believed that the chaotic behaviour generating 
from nonintegrability is essential criterion for realising 
Fourier's law in classical systems. Analogously from the 
numerical study of quantum systems with the coupling 
to external heat baths [l^, EH or without , it is ar- 
gued that the emergence of diffusive behaviour is related 
to onset of quantum chaos. Now we try to analyse the 
underlying mechanism of getting diffusive behaviour in 
this self-consistent reservoir model. Originally Bolsterli, 
Rich and Visscher [6| proposed the self-consistent reser- 
voir model to incorporate phenomcnologically the inter- 
actions of phonons with other degrees of freedom such 
as electron's charge and spin present in the physical sys- 



tem. Due to stochastic interactions with the internal 
reservoirs the inherently non-ergodic harmonic chain be- 
comes ergodic. Here the self-consistent reservoirs pro- 
vide the mechaism of scattering for phonons which is 
very much essential to get diffusive behaviour. It can 
also be posed in a different way that the self-consistent 
reservoirs act as the environment in a persistent random 
walk of phonon in a lane and break down the coherent 
nature of transport. In this context this model is sim- 
ilar to the models of particle transport studied in 11 1 
again with self-consistent particle reservoirs and in 15 1 
with heat baths modelling the dissipative environment. 
Now we point out certain inconsistencies in the applica- 
tion of "quasiclassical approximation" in Ref.0 which 
treats system classically neglecting all quantum fluctu- 
ations and random noises from the baths as quantum 
mechanically correlated. At high temperatures, where 
thermal fluctuations predominate over quantum fluctua- 
tions, the system is inherently classical. In the opposite 
limit, the strength of the anharmonicity in the quartic 
onsite model is weaker if the temperature is lower. Here 
the anharmonicity can be treated perturbatively in an 
effectively harmonic system. So in these two limits the 
so called "quasiclassical approximation" is valid. But 
in the intermediate temperature where the anharmonic- 
ity has significant strength, quantum fluctuations due to 
non-commutativity of the operators play much impor- 
tant role in lower dimensions. Then the "quasiclassical 
approximation" ceases to be correct. Thus, though use 
of the "quasiclassical approximation" looks attractive, it 
has probably little application for real problems of quan- 
tum transport where phonon-phonon interaction is cru- 
cial. Finally one main feature of our analysis is that the 
effective transport mean free path distinguishing ballistic 
regime from diffusive one, emerges naturally in the study. 
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